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ABSTRACT 

We have obtained 1087 near-infrared (NIR; JHKg) measurements of 21 Type 
la supernovae (SNe la) using PAIRITEL, the 1.3-m Peters Automated InfraRed 
Imaging TELescope at Mount Hopkins, Arizona. This new set of observations 
nearly doubles the number of well-sampled NIR SN la light curves. These data 
strengthen the evidence that SNe la are excellent standard candles in the NIR, 
even without correction for optical light-curve shape. We construct fiducial NIR 
templates for normal SNe la from our sample, excluding only the three known 
pecuUar SNe la: SN 2005bl, SN 2005hk, and SN 2005ke. The i/-band absolute 
magnitudes in this sample of 18 SNe la have an intrinsic RMS of only 0.15 mag 
with no correction for light-curve shape. We found a relationship between the 
i7-band extinction and optical color excess of Ah = 0.2E{B — V). This variation 
is as small as the scatter in distance modulus measurements currently used for 
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cosmology that are based on optical light curves after corrections for light-curve 
shape Combining the homogeneous PAIRITEL measurements with 23 SNe la 
from the literature, these 41 SNe la have standard if-band magnitudes with an 
RMS scatter of 0.16 mag. The good match of our sample with the literature sam- 
ple suggests there are few systematic problems with the photometry. We present 
a nearby NIR Hubble diagram that shows no correlation of the residuals from the 
Hubble line with light-curve properties. Future samples that account for optical 
and NIR light-curve shapes, absorption, spectroscopic variation, or host-galaxy 
properties may reveal effective ways to improve the use of SNe la as distance in- 
dicators. Since systematic errors due to dust absorption in optical bands remain 
the leading difficulty in the cosmological use of supernovae, the good behavior of 
SN la NIR light curves and their relative insensitivity to reddening make these 
objects attractive candidates for future cosmological work. 

Subject headings: distance scale - supernovae: general 



Introduction 



Type la supernovae (SNe la) observed at optical wavelengths in the rest frame have 
played a leading role in extragalactic astronomy and cosmology in the past decade. SNe la 



have been the ke y to measuring the Hubble constant (IFreedman et al 



2001 



Jha et al. 



19991: 



Riess et al. 2005 ) and demonstrating cosmic acceleration ( Riess et"al] 19981 : Perlmutter et al. 
I999I ). Thlse discoveries rest on a fou ndation of photometric and spectroscopic similarities 



between high- and low-redshift SNe la (jHamuv et al. 



2001; 


Hook et al. 


2005 


2007; 


Gare et al. 


2007 



31ondin et al. 



200 



Bronder et al 



i 



2007 



ha et al 



996 



■ liess et al. 



2006|: 



Foley et al. 



1999 



Conlev et al.ll200 



200 



a 



of cosmic expansion ha s been traced out to redshift z ^ 1 (ITonry et al.l l2003l : iKnop et al. 



Goldhaber et al 



Ellis et al. 



Garavini et al 



20071). The history 



2OO3I : iBarris et al.l |200J) . The predicted turnover to the matter-dominated era and the cor- 



respo nding deceleration at 2; > 1 have been seen by space-based work from HST ( iRiess et al. 



20071 ). Tw o current projects: ESSENCE (Equation of S tate: SupErNovae trace Cosmic 
Expansi on: iMiknaitis et al.l 120071 ; IWood-Vasey et al.l 120071 ). and SNLS (SuperNova Legacy 
Survey; lAstier et al.l l2006l ) seek to constrain the nature of the dark energy responsible for 
cosmic acceleration by measuring the equation-of-state parameter w with large, homoge- 
neous samples of SNe la. ESSENCE and SNLS data are both consistent with dark energy 
described by a cosmological constant, w = —1.0 to better than ±0.1. 



These cosmological results are based on measurements of optical emission in the super- 
nova rest frame, but recent work suggests that SNe la may be superior distance indicators 
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in the near infrared (NIR), with a i iarrow distribution of peak JHK,^ magnitud es and 4-8 
times less sensitivity to reddening (IMeikld l2000l : iKrisciunas et al.l l2004al . 120071 ). Redden- 
ing and absorption of supernova hght by dust in th eir host galaxies poses the single mos t 
vexing systematic question facing SN la cosmology (IWang et al.ll2006l : IConley et al.l 120071 ). 
Avoiding these complications would be very desirable, if the properties of the supernovae 
permit. These questions can be approached by building a homogeneous and large set of NIR 
observations. This paper is a step in that direction. 



Early photornetric NIR observations of S Ne la were made b v iKirshner et al.l (1l973l ): 



Elias et all fll98ll. Il985r): iGraham et all (119881): 



light curves publish e d bv iJha et al.l (I1999lll 



Frogel et al 



Fernandez et al. 



Candia et al. (12003): Krisciunas et al. (200G 



e: 



J2004bl . cl.l2005il2006h 
2006h: IPastorello et al.N2007bm: 



ias-Rosa et al.l (120061. 



2001 



Prieto et al 



200 



20031) 



1987). with more recent JHKg 



(2000): Valentini et al.l (2003) 



Benetti et al. ( 2004 ): 



Krisciunas et al. 



mm-. 



Krisciunas et al 



(120071) :IStritzinger k SoUermanl (120071 ) : IWang et al 



Phillips et al. 



: IPignata et al. ( 2008 ): Taubenberger et al. ( 20081 ). Northern hemisphere work from 



the CfA Supernova Program uses observa tions from the robo tic 1.3-m Peters Automated 



InfraRed Imaging TELescope (PAIRITEL; iBloom et al.ll2006bl ) at Mount Hopkins, Arizona 



The NIR observations are part of a system atic program of supernova observations that in 



eludes dense sampling of superno va spectra (lMathesonll2007l ) and concurrent UBVri optical 



photometry (iHicken et al.l l2008l . in prep.). This paper is a first report on the NIR data 



on SNe la from PAIRITEL. Similar work is underway in the southern hemisphere by the 



Freedman 


2005; 


Hamuv et al. 


2006) 



with progress in modeling SN la NIR light curves (IKasenl 120061 ) raises hopes that SNe la, 
especially in the rest-frame H band, can be developed into the most precise and accurate of 
cosmological distance probes. 



Recent work by iKrisciunas et al.l (j2004al ) shows that SNe la have a narrow range of lumi- 
nosity in JHKg a t the time of -B-band max imum light (t^max) with smaller scatter than in the 
B and V bands. iKrisciunas et al.l (l2004a] ) found no correlation between optical light-curve 
shape and intrinsic NIR luminosity. The NIR behavior is in sharp contrast with the optical 
light curves, where a variety of ingenio us methods have been devised to reduce t 
distance estimates: the Amis method ( Hamuy et al. 1996; Phillips et al. 19991: 



le scatter in 



20Q&\. the mult icolor light-curve shape method (MLCS/MLCS 2k2 



Jha et al.ii2007r ). the "stretch" method (iPerlmutter et al.iil997i: iColdhaber et al.li200lr) . the 



adap tive template method (SALT/SA LT2; iGuy et al.l i2005l : iAstier et al.l i2006i : iGuy et al. 



Riess et al 



'rieto et al. 



199 



a 



1998i : 



color-magnitude intercept calibration meth od (CMAGIC; Wang et al. 2(303h. the spectral 



2007i). and SiFTO (jCorilev et al.i 120081 ). For a sample of 16 SNe la observed in NIR bands, 
Krisciunas et all J2004ah found a RMS of aj = 0.14, an = 0.18, and ax, = 0.12 mag. We 
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set out to construct an independent set of observations to test this remarkable result. 



In this paper we present 21 SNe la observed with PAIRITEL from 2005-2007. The 
1087 individual data points in this sample represent the largest homogeneous set of SNe la 
to date and double the NIR SN la observations in the literature. Data collection with 
PAIRITEL is discussed in ^ Data analysis, including the mosaic creation process and 
our photometry pipeline, is discussed in §31 We construct NIR JHKg templates from this 
new sample as detailed in §1] and fit these templates to the SN la light curves to derive the 
JHKg magnitudes at the time of 5-band maximum light. These magnitudes are remarkably 
uniform, particularly in the J and H bands. Distances are calculated in §5], and we discuss the 



Krisciunas et al. 


(2C 


04i 


J ) as well as to some more recent SNe la ( 


Krisciunas et al. 


2005, 


2006, 


2007; 


Phillips et al. 


2006 


Pastore 


b et al. 


2007a; 


Stritzineer & SoUerman 


20071; 


Wane et al. 


2008; 


Elias-Rosa et al. 


20061 


2007 


Pignata et al. 


20081: 


Taubenbereer et al. 


2008|) in m\ and 



show that the J and H bands at maximum light, even from this heterogeneous sample, have 
a small range in absolute magnitude. Our conclusions are summarized in §81 



2. Data 



We present 1087 near-infrared (NIR) measurements with a signal-to-noise ratio > 3 of 
21 nearby SNe la obtained from 2005 to 2007 using the f/13.5 1.3-m PAIRITEL at the F.L. 
Whipple Observatory on Mount Hopkins, AZ. (Tabled]). Dedicated in Octob er 2004, PAIRI 



TEL uses the northern telescope of the Two Micron All Sky Survey (2MASS; ISkrutskie et al. 



20061 ) together with the 2MASS southern camera. PAIRITEL is a fully automated, robotic 
telescope with the sequence o f observations controlled by an intelligent queue-scheduling 
database (IBloom et al.ll2006bf). Two dichroics allow simultaneous observing in JHKs (1.2, 
1.6, and 2.2 nm: ICohen et al.l[2003[ ) with three 256x256 pixel HgCdTe NICM0S3 arrays. 
The image scale of 2"/pixel provides a field of view of 8.53' x 8.53' for each filter. Since the 
supernova observations are conducted with the instrume nt that defines th e 2MASS photo- 
metric system, we use the 2MASS point source catalog (ICutri et al.ll2003l ) to establish the 
photometric zeropoints. Images are obtained with standard double-correlated reads with 
the long (7.8 sec) minus short (51 msec) frames in each filter treated as the "raw" frame in 
the reduction pipeline. The telescope is dithered (< 2') every fourth exposure to aid with 
reductions. 

The CfA Supernova Program used PAIRITEL to follow up SNe discovered by optical 
searches at 5 > —30 degrees with ^ ^ 18 mag. These are objects for which we routinely 
obtain extensive data sets of optical spectra and CCD photometry using other instruments 
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at Mount Hopkins. With ~ 30% of the time on this dedicated robotic telescope avail- 
able for supernova observations, we observe ~ 3-5 SNe per night, and obtain host galaxy 
template reference images for each supernova. The first p ublished PAIRITEL s upernova 



observations were of the unusual core-collapse SN 2005bf (iTominaga et al.ll200"5[). There 



have 


)een previous papers that present the data processing from PAIRITEL 


Bloom et al. 


2006a 


Butler et al. 


2006 


Yost et al. 


2006; 


Bloom et al. 


2007; 


Blake et al. 


2007; 


Eisner et al. 


2007 




Modjaz et al. 


2008 


1 , but since NIR photometry remains more challenging than optical 



photometry we present a brief sketch of the data analysis path in the next section. 



3. Data Analysis 



Data were processed with one of two mosaicking pipelines. For fields dominated by large 
galaxies we used a pipeline (see lModjag|2007l : iKocevski et al.ll2007l ) that makes use of a bank 
of existing sky frames. For all other fields, we used the standard pipeline (see §3.11) . Both 
pipelines process the 100-150 images, each of 7.8 seconds duration, by flat correction, dark 
and sky subtraction, registration, and stacking to create a final mosaicked image in each of 
the JHKs filters. Mosaicked images typically comprise 3 image s at each dith er position, with 
a mosaicked FOV of ~12'xl2'. The raw images are SWarped (jBertinll2005l ) using 0.1"/pixel 
sub-sampling into final mosaics with l"/pixel sampling. Typical 1800-second observations 
(including slew overhead) reach magnitude limits of ~ 18, 17.5, 17 mag for J, if, and Kg 
respectively. 



3.1. Mosaics 

The PAIRITEL camera has no shutter, so dark current cannot be measured indepen- 
dently, and "skark" background frames include both sky and dark photons. PAIRITEL 
supernova observations did not include pointings that alternate between the source and the 
sky, so skark frames were created for each mosaic. 

For large host galaxies with angular size > 2' (in the 12' FOV), host galaxy contamina- 
tion prevents a reliable calculation of the skark background calculated from a pixel-by-pixel 
median through the stack of dithered images. Here skark frames are estimated using a me- 
dian match to an archive of skark frames in relatively empty fields observed on the same 
night. The archival skark image is selected as close in time to the supernova observations as 
possible, with a median matched in a 50x50 pixel box, known to have stable pixel properties, 
in the lower left quadrant of the array. After the dark plus sky frame is subtracted, images 
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are mosaicked using the drizzle technique (IFruchter fc Hookll2002l ) 



For well-isolated supernovae with host galaxies of angular size < 2'(in a 12' FOV), skark 
frames can be constructed from the science images themselves. This construction is done by 
applying a cubic B-spline spline fit to a complete time series of sky values in each pixel for 
that observation. Since the fitting characteristics of the spline curve vary by filter, elevation, 
and weather, parameters are adjusted to best fit each observation and filter. The pixels 
with sources are masked out using object masks generated from Source Extractor catalogs 
extracted from the raw frames. The pipeline replaces each masked pixel's flux value with 
a more simple "median skark" pixel value that is the sum of the median flux value of each 
frame and a median over time of the pixel's deviation from this median frame value. 

The spline fit parameters define a function that tracks the time variation of the skark 
background in each of the 256x256 pixels over the duration of the dither pattern. We 
construct a skark image at the time of each 7.8-second science exposure by calculating the 
skark value in each pixel at this time and subtract that skark image from the science image. 
For cases where the galaxy occupies a small fraction of the field, this process improves the 
sky subtraction. This standard pipeline works better than simply taking the pixel-by-pixel 
median through the dithered (unregistered) stack of science images because it samples the sky 
variation on time scales th at are shorte r than the total exposure time. Mosaicked images are 



constructed using SWarp (lBertinll2005l ). We use SWarp with a simple BILINEAR sampling 



because such a kernel is most suited to our under-sampled raw images. 

Most of our objects were reduced with the standard processing code. Only SN 2005cf, 
SN 2005ke, SN 2005na, and SN 2006X required the large galaxy approach. For both methods, 
bad pixel masks and fiat fields were created from archival images. Fig. [1] shows a final JHKg 
color mosaic created in the standard processing mode. 



3.2. Photometry 



Mosaicked images were fed to the photometry pipeline we have used in the ESS ENCE 
and SuperMACHO projects jRest et allbood : ICarg et al.lboOTi : iMiknaitis et al.lboOTh . This 
pipeline determines the photon noise for each pixel based on the sky noise in the mosaic im- 
ages, registers the images for a supernova to a commo n reference frame, and performs point- 
spread f unction (PSF) ph otometry using DoPHOT (ISchechter et al.lll993l ). The 2MASS 
catalog (ICutri et al.ll2003l ) is the natural astrometric and photometric reference system for 
these observations. In a typical 12'xl2' field of view, there were 10-100 2MASS stars in each 
filter. These stars were sufficient to calibrate the images to the 2MASS system JHKg firmly 
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enough that the underlying uncertainty in the 2MASS system, of about 3%, is the domi- 



nant e rror in the photometric cahbrations for the hght curves presented her e . iCohen et al. 
fl2003h describe the 2MASS JHR^ fiher system in detail and Iheggett et all (|2006|) provide 



color transformations to put observations with the 2MASS filter system on other widely-used 
photometric systems. 

Final reference images were taken for each galaxy after the supernova had faded. DoPHOT 
photometry on the images was used for those SNe la that were clearly separated from their 
host galaxy and had little underlying contaminating light (SN 2005ao, SN 2005cf, SN 2005el, 
SN 2005hk, SN 2005ke, SN 2005eu, SN 2005iq, SN 2005na, SN 2006N, and SN 2006X). The 
standard background annulus subtraction used by DoPHOT to estimate the local sky was 
sufficient to remove any remaining contaminating galaxy light. In every case, the contribu- 
tion from underlying galaxy light was less than 10% of the SN la light at maximum. Seeing 
at PAIRITEL is limited by the dome seeing and remains relatively constant from 1.8-2.0". 
The good effect of this mediocre dome seeing is that the contaminating galaxy light within 
the PSF is very nearly constant. 



We used subtraction-based photometry (following iMiknaitis et al.l (120071 ) ) for SNe la 



that w ere not clearly separate d from their host gal axy. The NN2 method of iBarris et al. 



(120051 ) was employed (as used in IMiknaitis et al.ll2007l ) by subtracting all Nx[N~l)/2 unique 
pairs of images to minimize the sensitivity to subtraction errors. Because our PAIRITEL 
data are not critically sam pled, the reliabili t y of t he PAIRITEL image subtraction is not as 
good as the subtractions in IMiknaitis et al.l (120071 ). Our photometric pipeline automatically 
detected and screened out subtractions with significant residual fiux in known stars, leaving 
us with fewer high-quality light-curve points than would have been obtained in the case of 
perfect subtractions. For the supernovae where the underlying galaxy contribution was a 
small fraction of the supernova light, direct photometry on the unsubtracted images proved 
preferable. 

Photometry was extracted from either the unsubtracted or the subtracted images by 
forcing DoPHOT to measure the PSF-weighted fiux of the object at a fixed position. This 
position was determined in the J-band difference images that had a signal-to-noise ratio > 5. 
This average position was used for fixed-position DoPHOT photometry of each image of a 
SN la. 

Flux measurements were calibrated t o the 2MASS system by using photometric solutions 
to the 2MASS catalogs (ICutri et al.ll2003l ). For the difference images the calibrated zeropoint 
from the template w as used, with su i table correction for the convolution of the template 
image as detailed bv IMiknaitis et all J2007h . Our NIR PAIRITEL light curves for the 21 
SNe la presented in this paper are given in Table [51 
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4. NIR Template 



We constructed JHKg templates with NIR light curves from PAIRITEL, using our own 
optical CCD observations (IHicken et al.l l2008l ) to establish the time of B-band maximum, 

In constructing the templates, we excluded the three fast-declining SNe la from our 
overall sample of 21 SN e la: SN 2005bl, SN 2005 hk, and SN 2005ke. SN 2005bl has an 
optical fast decline rate (ITaubenberger et al.ll2008l ) similar to the archetypal fast-decliner, 
SN 19 91bg JFihppenko et aljll992h . SN 2005hk is known to be an unusual SN la JphiUips et al 



2007; Sahu et al 



20o'7fi. SN 2005ke is a 199 1bg-like SN la with possible circumstellar inter- 



action flPatat et al.ll2005l : llmmler et al.ll2006h . SN 2005bl, SN 2005hk, and SN 2005ke exhibit 
only one infrared hump instead of the two exhibited by normal SNe la. The JHKg maximum 
for these "dromedary" light curves occurs after 5-band maximum light, unlike normal "bac- 
trian" SNe la for which the first NIR maximum occurs 3-5 rest-frame days before 5-band 
maximum light. We compare our SN la if-band template to SN 2005bl, SN 2005ke, and 
SN 2005hk in Fig. HI Because we have excluded these three objects, our template does not 
extend to 1991bg-like SNe la or other unusual SNe la. This seems a reasonable approach 
while the data set is small. Unusual supernovae can be identified from a NIR light curve 
alone, or they can be identified from their optical spectra and light curves which we obtain 
as a matter of course. In any case, we defer incorporating these unusual objects into a more 
comprehensive treatment of NIR light curves until the available database of near-infrared 
SN la light curves is more fully populated. 



We apphed the i^'s-corrections of iKrisciunas et al.l (l2004bl ) to the NIR hght curves based 
on the heliocentric redshift. The NIR light curves for the remaining 18 SNe la were reg- 
istered to a common phase by su btracting ternax and accounting for time dilati on based 
on the heliocent r ic reds hift, ^heiio ( jLeibundgut et al.lll996l : iGoldhaber et al.ll200ll ). Unlike 
Krisciunas et al.l (l2004al ). we do not further adjust this phase by the optical light-curve width 
parameter. We found that compensating for the light-curve width gave no improvement in 
the H-band template fits (which is the focus of this paper) or in the resulting dispersion of 
absolute magnitudes (see ^ determined between —10 and -|-20 days in phase. However, 
the position of the J-band second NIR maximum is variable, and may be related to intrinsic 
luminosity. In larger data sets and in more detailed investigations of the J band it will 
be worth exploring whether a width-correction parameter would produce a more effective 
template for fitting the second maximum at later times. 

Between [—10, -|-30] days, where we have excellent sampling, we constructed our lightcurve 
templates by an unweighted averaging of the points from this initial set of standardized light 
curves into 1-day bins. No outlier rejection was used in this process. The SN la light curves 
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were then registered to a common magnitude scale by fitting each to this zeroth-iteration 
template and determining the maximum magnitude at B-band maximum light (IRemax)- 
The fitting process is insensitive to the determination of the initial IRsmax because this pro- 
cess is repeated with each iteration and the memory of the first guess is completely erased 
after the second iteration. At late times (> 30 days), we had fewer points from which to 
construct the template, but it was clear that the SNe la were consistent with a linear fit. 
For our templates after +30 days we adopted hnear templates based on a fit to SN la data 
at these late times with a forced connection to the binned value at 30 days. The complete 
template was then smoothed twice with a boxcar of length 5 days. The smoothing has very 
little effect near maximum light because of our dense sampling at those epochs. Varying the 
smoothing length between 1 and 5 days did not substantially affect the template. 

To refine this template, we determined the magnitude value at ismax that minimized the 
of the template versus the data over the span from [—10, +20] days. These new magnitude 
offsets were then used to seed the process above where we had previously used zero for the 
IR-Bmax values. The procedure was iterated twice more to construct final values at t^max for 
each SN la, and the final template was constructed based on these values. We found that 
three iterations were sufficient to reach convergence. The uncertainty in the template is the 
standard deviation of the residuals of the SN la light-curve points around the mean template 
in a moving 5-day window. 

We also explored using a simple linear interpolation on 1-day sampling and found that 
the binned template and the linearly-interpolated templates agreed to within hundredths 
of a magnitude in this region near peak. We take this as a strong indication that we have 
reached limit of intrinsic variabihty of SNe la in the construction of this template since the 
effective use of 5 times more data in the binned approach yielded the same template as the 
hnear interpolation, which only uses ~ 20% of the available data. 

This template construction procedure was performed for each of the JHKg NIR pass- 
bands. For our data, the if-band light curve was observed to exhibit small scatter from 
[—10, +30] days with a particularly tight distribution from [—10, +20] days. The J-band 
aggregate hght curve has small scatter about the template from [—10, +10] days but begins 
to show variations in the time and flux of the secondary maximum for individual supernovae. 
The Kg-haxid data from our PAIRITEL observations is not as good as the J- and i^-band 
data due to increased sky background in Kg. Our template in the Kg band shows signif- 
icantly more scatter with fewer objects. We draw no firm conclusions in this paper with 
regard to i^^-band SN la light curves. 

The unusual cases that were excluded from the template construction are easily seen 
to be far from the H-band template. Fitting SN 2005bl, SN 2005hk, or SN 2005ke with the 
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PAIRITEL SN la iJ-band template results in a sufficiently bad x^/DoF (> 3) (see Table [3]) 
that these SNe la would be excluded based on that criterion alone. 



5. Galaxy Redshifts 

Recession velocities for our SNe la are those of the identified host galaxies in the 
NASA/IPAC Extrag alactic Database (NE pfl We use the quoted velocity with respect to the 
Virgo infall model of iMould et al.l (120001 ) . For SNe la within 3000 km s^^ we used available 
information on distances to t he galaxies to minimize o ur sensitivity to the details of the local 
flows. We u s ed th e values of iKrisciunas et al.l J2004al jbl) for SN 1998bu and SN 2002bo and 
Wang et al.l (120081 ) for SN 2006X. After accounting for these SNe la, we were left with three 
SNe la within 3000 km s^^ for which we used the IMouM et all J200oh values, SN 2003du, 
SN 2005cf, and SN 2006D. Outside of 3000 km s~^ differences among various flow mod- 
els are not significant. One supernova, SN 1999cl, occurred in a Virgo cluster member 
galaxy, NGC 4501, so we used the Virgo cluster redshift instead of the recession velocity 
of NG C 4501. We assumed a peculiar velocity of a^ei = 150 km s~^ (IRadburn-Smith et al. 
20041 ) and the redshift measurement uncertainty for each host galaxy as reported by NED. 
These uncertainties were converted into a distance modulus uncertainty, (j^(vel), using Eq. [H 



0-, 



(vel) 



zlnlO 



This is valid for a smooth cosmology that resembles any currently allowed cosmological 
expansion model. Values for cr^(vel) for the SNe la considered in this paper are given in 
column 12 of Table O 

We confirmed the NED heliocentric redshift values with the heliocentric galaxy redshifts 
derived from spectra of the hosts as obse rved by the CfA Supe rnova Program^ and measured 
using the galaxy redshift code RVSAO (iKurtz &: Minklll998l ). All were consistent with the 
NED values within 200 km s~^. For SN 2007cq, no NED redshift was available for the host 
galaxy so we used the value frorn the CfA galaxy spectrum as the heliocentric velocity and 
then used the Mould et al.l (120001 ) model conversion provided by NED to express this redshift 
in the CMB+Virgo infall frame. 



http : //nedwww. ipac . caltech. edu/j 



' http : //www . cf a . harvard . edu/ supernova/RecentSN . html 
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6. Distance Moduli and the Standard NIR Brightness of SNe la 



The most uniform behavior of our SN la hght curves is in the if-band, and that band 
is where we focus our analysis. In the most simplistic treatment, we take the -ffBmax for each 
SN la to have a standard value. Then the distance modulus is simply: 



Bmax 



- M 



(2) 



where M^^^^^ is the absolute if-band magnitude of a SN la at i?-band maximum light. The 
numerical value for the absolute magnitude of SNe la in H is separate from the question of 
how well the if-band magnitudes trace cosmic expansion. 

When we combined the results for the PAIRITEL SNe la with our analysis of the litera- 
ture SNe la in §7] we found that an if-band extinction correction of = 0.19E{B — V) best 
minimized the correlation between Hubble-diagram residual and E{B — V) as measured by 
MLCS2k2 and applied that correction to the apparent magnitudes computed from fitting to 
our if-band template. While, the trend was not formally statistically significant because the 
most highly extinguished SNe la are also the ones that are at the lowest recession velocity 
and have the most uncertain distances it was consistent with direct comparison of extinction 
for SNe la in the same galaxy. Two of the SNe la in our sample, SN 2002bo and SN 2002cv, 
were both in NGC 3190 and had significant extinction. When we assumed that each of these 
SNe la had the same intrinsic brightness, we found the same relationship between if-band 
extinction and optical color excess of Ah = 0.2*E{B — V). While quantitatively uncertainty. 



this Rh ~ 0.2 corresponds to an equivalent Ry = 1.2, which is in qua^ 
previous work on observed supernova color- luminosity behavior (c.f.. 



itative agreement with 



Perlmutter et al. 1999 



Knop et aklboosl : IGuv et al.lEoOsI : IConlev et allboO?! : IGuv et al.l[2007h . Resolving this mys- 



terious behavior of supernova colors, host galaxy dust, and apparent luminosity is currently 
the major challenge in SN la cosmology. 

The set of redshifts, apparent magnitudes, and accompanying uncertainties from ^were 
compared to a ACDM concordance cosmology (^m, ^a, w) = (0.27, 0.73, — 1). We solve for 
the absolute mean brightness of a SN la in H-band at B-band maximum light by taking 
the weighted mean of the differences of the apparent magnitudes and the ACDM distance 
moduli (the exact cosmology doesn't quantitatively affect the results at these redshifts): 



N 



TTi . i 

-"Bmax P 




1 



(3) 



Assuming an Ho of 72 km s ^ Mpc ^ (IFreedman et al.ll200ll : ISpergel et al.ll2003l . 120071 ) we 
find Mur, = —18.08 ±0.04 with an RMS of 0.15 mag for the fiducial if-band magnitude at 

J-* Bmax O o 
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B-band maximum light. A similar analysis for the J-band magnitude at 5-band maximum 
light finds Mjg^^^ = —18.29 ± 0.09 with an RMS of 0.33 mag. These values are for the new 
SNe la presented here. 

For a cosmological fit one would marginalize over the combination 

-^//B.a. = Mj^3_-51ogio/io + 15 + 51ogioc (4) 
= S^logioe (5) 

(6) 

where h = Ho/{100 km s~^ Mpc~^) and c is the speed of light in vacuum in km s~^. The 
marginalization in our case is trivial because we are assuming a fixed ACDM cosmology. 
This assumption allows us to propagate a Ho uncertainty of 10% directly through Eq. [6] to 
determine our global uncertainty in Mhj^^^- The mean value of Mh^^^^ is unchanged by 
the marginalization, but the total uncertainty is increased to ±0.22. The uncertainty in Ho 
dominates the overall uncertainty of M^jj^^^^^^- 

Small variatio n in the NIR absolute ma gnitudes of SNe la is an important result that 



was first found by iKrisciunas et al.l (j2004al ). Here, we have an independent confirmation 
of this result from a homogeneous set of well-sampled SNe la light curves observed with 
PAIRITEL. In the next section, we expand that data set by joining the PAIRITEL sample 
with SNe la from the literature. 



7. Comparison with Literature SNe la 



We joined our PAIRITEL sample with 23 SNe la from the literature (see Table [3]) ex- 
tending back to SN 1998bu (earlier observations of supernovae with earlier detectors were 
not included in this study). When creating the template including t he literature SNe la , 
SN 1999ac was excluded because it was classified as an unusual SN la (jPhillips et al.ll2006l ). 
The PAIRITEL SNe la listed in § H] were similarly exclude d from the global te mplate- 
generation. The literature SNe la that were published on the Persson et al.l (1998 ) system 
were converted to the 2MASS system using the conversion equations of ICarpenterl (120011 ). 
These corrections were typically on the order of 0.01 mag. We generated a new set of NIR 
templates following the same prescription as in § H] and found that the full sample of SNe la 
was just as well-behaved as the PAIRITEL sample. This is a tribute to the care of the ob- 
servers and the clear definition of the 2MASS photometric system. The NIR SN la templates 
and associated standard deviations for this combined sample are shown in Fig. [3] and given in 
Table [21 The small scatter of the H-hand light curves from [—10, +20] days and the J-band 
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light curves from [—10, +10] days was preserved in the heterogeneous sample. Fig. [5] shows 
the tight scatter of the if-band SN la light-curv e data points around the fiducial template. 
Over- plotted in Fig. [3] are the NIR templates of iKrisciunas et al.l (l2004al ). There is general 
agreement but noticeable systematic discrepancy in the H-band template before maximum 
light. Together with the difference in definition of quo ted if -band max i mum ( in this pa- 
per we refer to Mh at B-band maximum light whereas IKrisciunas et al.l (l2004al ) quote the 
observed if-band magnitude at H-band maximum light), this disagreement in the if-band 
templates explains the bulk of the 0.3 mag difference between the quoted Mh values in these 
respective works. 

We present the if-band apparent magnitudes from this global template in Table |3l The 
comparison of these apparent magnitudes with the expected luminosity distance relation of 
a ACDM cosmology, (^m, ^^a, K,^) = (0.23,0.77,0.72,-1), is shown in Fig. El We find 
a global if-band RMS of 0.16 mag around an absolute magnitude of Mh = —18.07 ± 0.03 
(see Fig. [7|). A similar analysis for J-band finds similar agreement with the PAIRITEL-only 
values with a global value of Mj = -18.27 ± 0.07 with an RMS of 0.29 mag. The i^^-band 
data from PAIRITEL is more limited and only adds 9 normal SNe la to the literature. The 
updated global values including these new SNe la against the global template as generated 
for J and H are Mk = -18.30 ± 0.11 with an RMS of 0.29 mag. To check the quality 
of the light-curve fits and to see if the fits are robust, we plot the H-hernd residuals with 
respect to the x^/DoF and degrees of freedom of the fits to the if-band template in Fig. [HI 
Fig. [H] demonstrates that after an extinction correction of Ah = 0.19E{B — V), there is no 
correlation in the ii-band residuals with the light-curve sha pe parameter, A , or extinction. 
Ay, as measured from the optical light curves by MLCS2k2 (jjha et al.ll2007l ). 

To reduce the effects of uncertainties in the local flow, we also calculated the RMS 
for only the 28 SNe la with a CMB+Virgo infall model velo city of > 2000 km s~^ We 
found an RMS of 0.15 mag for this sample. Previous work by iRadburn-Smith et al.l (120041 ) 



and iNeill et al.l (120071 ) found that SN la distances were more consistent with the Hubble 



law when a local flow model was incorporated. Future work from larger samples of nearby 
SNe la in the NIR will allow for the investigation of both the local flow and the absolute 
NIR magnitudes of SNe la from integrated analyses with more sophisticated flow models of 
the local Universe. 



8. Conclusions 



We have constructed improved NIR templates to find accurate luminosity distances 
to SNe la using the if-band light curve from [—10, +20] days after 5-band maximum light. 
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Within the photometric and local-flow uncertainties, the distribution about a best-fit ACDM 
model is consistent with no intrinsic dispersion of our if-band template fit magnitude values. 
Even without accounting for peculiar velocity uncertainties, the RMS scatter of 0.15 mag 
from these uncorrected inferred luminosity distance moduli is as small as corrected distances 
from optical-based methods. No correction for light-curve shape has been made here and 
the treatment of dust extinction is relatively limited in its extrapolation over a factor of 2 in 
wavelength. It is quite plausible that using information on the NIR and optical light-curve 
shapes, dust absorption as fit from the optical to the infrared, spectroscopic variation, or 
host-galaxy properties will result in even smaller scatter and better distance determinations. 



The iKrisciunas et al.l (l2004al ) result was determined from an inhomogeneous sample of 
16 light curves. This paper improves on that work by providing a homogeneous sample 
of 18 PAIRITEL SN la light curves that more than doubles the number of independent 
photometric observations in the literature. With our ongoing program at Mt. Hopkins adding 
~ 10-15 SNe la per year, PAIRITEL will continue building an extensive nearby training set 
of ground-based NIR SNe la. While we observe all the SNe la that meet our criteria, objects 
in the Hubble fiow observed with adequate integration times to beat down the measurement 
errors will be the most helpful. 

In this paper, we considered reddening values derived from optical data, and found a 
relationship of apparent Ah = O.OGAv between the absolute H-hand magnit ude and the opti- 



cal re ddening. However, recent work on the highly extinguished SN 2002cv (lElias-Rosa et al. 



20071 ) shows that reddening corrections will be necessary for some highly-absorbed SNe la 
in the NIR. Reddening uncertainties derived from optical data alone currently represent the 
most signi ficant systematic error affecting SN la luminosity distance measurements. Simu- 
lations by IKrisciunas et al.l (120071') dern o nstra te major improvements to be gained from the 



addition of NIR data. IKrisciunas et al.l (120071 ) show that simulated distance modulus errors 
are improved by factors of 2.7 and 3.5 by adding J, and JHR^ to UBVRI data, respectively. 
Reducing systematics due to reddening are crucial to future space-based SN la surveys which 
will be large enough to avoid limitations from sample size. In that case, improving the con- 
straints on cosmological parameters will be limited by systematics. We will explore the 
utility of a full c ombination of optical and NIR data to more precisely measure and correct 
for reddening in iFriedman et al.l (120081 ) . 



Although ground-based NIR data can be obtained for low redshift objects, rest-frame 
NIR observations for high-redshift SNe la will have to be done from space. Currently, rest- 
frame S N la Hubble diagrams of high-redsh ift SNe la have yet to be constructed beyond the 
/ band (iFreedmanI l2005l : iNobili et al.l |2005|) , with l imited studies of SNe la and their hos t 
galaxies conducted in the mid- infrared with Spitzer (IChary et al.ll2005l : iGerardy et al.ll2007l ). 
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Because nearby SNe la are excellent standard candles in the NIR, it may be worth careful 
consideration of the rest-franae NIR for the space missions u sing SNe la for cosmology. For 
example, the SNAP jAlderinelbooi ) and DESTINY satellites JLauerll2005[ ). candidates for the 
NASA/DOE Joint Dark Energy Mission (JDEM) mission, are both currently designed with 
detectors sensitive out to 1.7 /im, which will only detect rest-frame if-band light (1.6 /im) 
out to z ~ 0.1. Only a detector capable of observing rest-frame if-band at z ~ 0.5-2 could 
take full advantage of Nature's gift to us of a superb standard candle in the rest frame H- 
band. Such a detector would require sensitivity from 2-5 /im. We can look forward to work 
of this type with JWST, just over the horizon. The astronomical community should discuss 
which aspects of a JDEM mission could best be carried out in the rest-frame NIR. 
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Fig. 1. — PAIRITEL JHKg composite color image of SN 2006D at a week past maximum 
light when the SN la had NIR magnitudes of {J,H,K) ~ (16.1,14.8,15.0). The image 
shown is 7. 5' x 7.5' in size. There are ~ 30 2MASS stars in this field that were used for the 
photometric calibration. The SI, S2, and S3 labels indicate three representative 2MASS 
stars with JHKs magnitudes of SI: (14.11, 13.76, 13.71), S2: (16.02, 15.335, 14.830), and 
S3 (13.23, 12.61, 12.41). 
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Fig. 2. — PAIRITEL NIR lightcurves of the 21 SNe la presented in this paper. There 
are 18 normal SNe la plus 3 fast-declining or unusual SNe la: SN 2005bl, SN 2005hk, and 
SN 2005ke. These SNe la are displayed on their own row at the bottom of the plot. 
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Fig. 3. — SN la JHKg templates. The 1-a uncertainties (dashed hnes) are based on the 
sample variance within a 5-day moving window of the SNe la with respect to the fiducial 
template. The vertical line marks the time after which a linea r fit was used to f it the data 
(0.3 mag uncertainties). The 3rd-order polynomial templates of iKrisciunas et al.l (j2004al ) are 
plotted in dash ed lines. The righ t -hand figure shows the difference between our template and 
the template of lKrisciunas et al.l (l2004al ). While the polynomial templates show approximate 
agreement with the new templates presented here, there is a noticeable difference in the H- 
band lightcurve before B-band maximum light. This di fference is largely respon sible for the 
different values of H-band maximum in this paper and IKrisciunas et al.l (l2004al ) . 
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T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 

H-band template 

Std. Deviation 




SN 2005bl + " 
SN 2005hk o 
SN 2005ke □ - 



^ I I I I I I I I I I I I I I I I I I I \ I I I I I I \ I I I 

-10 -5 5 10 15 20 

Restframe Phase [(Days - TQ^J/{^+z^^J] 



Fig. 4.— SN 2005bl, SN 2005hk, and SN 2005ke were observed as part of the PAIRITEL 
campaign but were excluded from the construction of the template because they were known 
to be unusual SNe la. The H-band template we have constructed for this paper is valid from 
-10 to +20 rest-frame days from B-band maximum light. We compare the if-band light 
curves of these unusual SNe la with the normal SN la if-band template to demonstrate the 
clarity with which these unusual supernovae can be distinguished from the normal SNe la 
used in the Hubble diagram. 
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T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — \ — I — I — I — I — I — I — \ — I — I — r 

H-band template 

Std. Deviation 




PAIRITELSNela • 
literature SNela o 

2 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

-10 -5 5 10 15 20 

Restframe Phase [(Days - TB^ax)/(1+Zheiio)] 

Fig. 5. — NIR if-band template based on all of the SNe la considered in this paper (excluding 
the dromedary SNe la shown in Fig. Hj) together with the constituent SN la observations. 
Model uncertainties are indicated by dashed lines as in Fig. [31 For clarity data error bars are 
not displayed. Instead points are scaled by their signal-to-noise ratio: small - 10 >SNR> 3; 
medium - 20 >SNR> 10; large - SNR> 20. 
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Fig. 6. — if-band SN la Hubble Diagram {fiH vs. z and residual vs. z). RMS dispersion 
of 0.15 mag. SNe la from the literature are shown as open diamonds while the new SNe la 
from this paper are shown as filled circles. The error bars are the fit uncertainties from the 
fit to the if-band template (column 7 of Table [3]). The reduced Hubble constant, ho, is 
K = Ho/ WO km s-i Mpc-^ 
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Fig. 7. — A histogram of the derived absolute magnitudes Mh^^^^ based on a ACDM cos- 
mology with Ho = 72 km Mpc~^for all SNe la (top panel) along with histograms for the 
PAIRITEL sample (middle panel) and the literature SNe la (bottom panel). Overlying each 
histogram is a Gaussian with the RMS width of the given sample (dotted line) normalized 
to the number of SNe la in each plot. 
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Fig. 8. — Hubble diagram residuals, (top panel) The Hubble diagram residuals of the 
extinction-corrected m^^ as a function of the light-curve shape parameter, A, and the mea- 
sured optical extinction, Ay- Both plots are are entirely consistent with no dependence of 
Hubble diagram residual on A after the extinction correction of Ah = 0.06^4^. (bottom 
panel) The Hubble diagram residuals as a function of x^/DoF and the number of H-band 
light-curve points (DoF=# light-curve points— 1). These plots demonstrate that the resid- 
uals are not sensitive to the formal x^/DoF of the fit and are reliable even with just a few 
data points. The large number of points with x^/DoF< 1 implies that we have likely overes- 
timated the uncertainties in our fiducial template or in our photometric observations. Unlike 
in Fig. El the error bars shown here represent both the model uncertainties (Table [31 column 
7) and the peculiar velocity (150 km s~^) (Table [31 column 4) added in quadrature (Table [31 
column 12) to better allow for a relative comparison of the deviation. There is a mild corre- 

Iflt.inn bfitwppn thp siyp nf thp" prrnr bars and At^ bfipaiise' thp SNp Ta with larcrp Ai/ alsn tp-nd 
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Table 1. SN la JHR^ Peak Magnitudes 





oN 




«^peak 




w b 

.'^peak 




peak 


M \T AT a 


Reference 


SN 


1998bu 


11 


55 


(0.03) 


11 


59 


(0.03) 


11 


42 


(0.03) 


18,20,20 


J99,H00 


SN 


1999cl 


12 


80 


(0.02) 


12 


77 


(0.04) 


12 


58 


(0.02) 


6, 5, 6 


KOO 


SN 


1999cp 


14 


50 


(0.02) 


14 


77 


(0.02) 


14 


57 


(0.06) 


2, 2, 2 


KOO 


SN 


1999ee 


14 


77 


(0.03) 


15 


00 


(0.03) 








33,33, 


K04a 


SN 


1999ek 


16 


10 


(0.05) 


16 


08 


(0.04) 


16 


54 


(0.07) 


14,15, 1 


K04b 


SN 


1999gp 


17 


31 


(0.06) 


17 


16 


(0.10) 


16 


51 


(0.14) 


3, 3, 3 


KOI 


SN 


2000E 


13 


47 


(0.02) 


13 


74 


(0.06) 


13 


42 


(0.06) 


16,15,10 


V03 


SN 


2000bh 


16 


80 


(0.03) 


16 


75 


(0.03) 


16 


64 


(0.04) 


20,21,15 


K04a 


SN 


2000bk 


17 


64 


(0.02) 


17 


45 


(0.05) 








18,16, 


KOI 


SN 


2000ca 


16 


36 


(0.03) 


16 


71 


(0.08) 








12,11, 


K04a 


SN 


2000ce 


16 


82 


(0.03) 


16 


12 


(0.03) 


15 


98 


(0.05) 


3, 5, 5 


KOI 


SN 


2001ba 


16 


99 


(0.02) 


17 


24 


(0.04) 


17 


06 


(0.06) 


15,14, 9 


K04a 


SN 


2001bt 


15 


34 


(0.07) 


15 


64 


(0.02) 


15 


40 


(0.03) 


21,21,11 


K04b 


SN 


2001cn 


15 


90 


(0.03) 


15 


66 


(0.05) 


15 


67 


(0.09) 


19,19, 3 


K04b 


SN 


200 Icz 


15 


40 


(0.06) 


15 


65 


(0.15) 


15 


40 


(0.12) 


12,12, 9 


K04b 


SN 


200 lei 


12 


93 


(0.02) 


12 


97 


(0.03) 


12 


84 


(0.05) 


33,32,32 


K03 


SN 


2002bo 


13 


68 


(0.02) 


13 


82 


(0.02) 


13 


88 


(0.04) 


14,14,14 


K04b 


SN 


2002CV 


14 


75 


(0.03) 


14 


34 


(0.01) 


13 


95 


(0.03) 


17,10,16 


DP02,ER07 


SN 


2002dj 


14 


56 


(0.09) 


14 


79 


(0.03) 


14 


50 


(0.09) 


21,21,21 


P08 


SN 


2003du 


14 


42 


(0.02) 


14 


66 


(0.02) 


14 


38 


(0.02) 


6, 6, 6 


St07 


SN 


2003cg 


13 


58 


(0.08) 


13 


68 


(0.10) 


13 


43 


(0.01) 


10,10,10 


ER06 


SN 


2004S 


14 


95 


(0.02) 


15 


04 


(0.02) 


14 


87 


(0.06) 


16,17,12 


K07 


SN 


2004eo 


15 


70 


(0.07) 


15 


79 


(0.07) 


15 


67 


(0.09) 


11, 9,10 


Pa07b 


SN 


2005ao 


17 


90 


(0.08) 


17 


65 


(0.00) 


16 


82 


(0.32) 


6, 6, 5 


WV08 


SN 


2005bl'' 


17 


35 


(0.04) 


17 


33 


(0.07) 


17 


47 


(0.14) 


15,12, 4 


T08,WV08 


SN 


2005cf 


13 


83 


(0.04) 


13 


96 


(0.02) 


13 


97 


(0.02) 


17,17,15 


WV08<= 


SN 


2005ch 


17 


48 


(0.06) 


16 


99 


(0.05) 


16 


20 


(0.14) 


10,10, 8 


WV08 


SN 


2005el 


15 


45 


(0.02) 


15 


55 


(0.03) 


15 


31 


(0.03) 


29,29,29 


WV08 


SN 


2005eq 


16 


74 


(0.03) 


17 


05 


(0.08) 


16 


51 


(0.11) 


30,26,27 


WV08 


SN 


2005eu 


17 


03 


(0.09) 


17 


20 


(0.16) 


16 


64 


(0.22) 


22,22,18 


WV08 


SN 


2005hkd 


15 


81 


(0.03) 


15 


60 


(0.05) 


15 


48 


(0.05) 


18,12,13 


WV08 


SN 


2005iq 


17 


30 


(0.09) 


17 


61 


(0.19) 








10, 3, 


WV08 


SN 


2005kc<^ 


14 


18 


(0.04) 


13 


98 


(0.03) 


14 


03 


(0.02) 


19,18,10 


WV08 


SN 


2005na 


17 


04 


(0.05) 


16 


96 


(0.14) 


16 


76 


(0.27) 


24,19, 7 


WV08 


SN 


2006D 


14 


32 


(0.02) 


14 


47 


(0.03) 


14 


54 


(0.05) 


22,22,20 


WV08 


SN 


2006N 


16 


16 


(0.11) 


15 


90 


(0.26) 


15 


72 


(0.32) 


14,12,11 


WV08 


SN 


2006X 


12 


86 


(0.03) 


12 


92 


(0.02) 


12 


74 


(0.04) 


60,66,47 


WV08'' 


SN 


2006ac 


16 


62 


(0.14) 


16 


60 


(0.20) 


16 


44 


(0.18) 


20,22,15 


WV08 


SN 


2006ax 


15 


77 


(0.05) 


16 


11 


(0.10) 


15 


89 


(0.06) 


17,16,16 


WV08 


SN 


2006cp 


17 


02 


(0.04) 


16 


96 


(0.15) 


16 


84 


(0.13) 


5, 5, 3 


WV08 


SN 


2006gr 


17 


84 


(0.07) 


17 


82 


(0.18) 


17 


72 


(0.25) 


6, 4, 1 


WV08 


SN 


20061e 


16 


15 


(0.10) 


16 


38 


(0.05) 


16 


05 


(0.17) 


38,32,36 


WV08 


SN 


20061f 


15 


59 


(0.08) 


15 


50 


(0.13) 


15 


24 


(0.14) 


30,33,26 


WV08 


SN 


2007cq 


16 


37 


(0.05) 


16 


92 


(0.21) 








5, 4,0 


WV08 



''Number of epochs w/ SNR> 3 in the JHKsVi^t curves, respectively. 

''Magnitudes of the maximum observed magnitude in JHKg. Magnitude errors are 1-a 
symmetric errors. 
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Re ference codes: WVOS (P AIRITE L photometry; thi s pape r); J99: Ijha et alj lll99g^; 



HOO: 



iHernandez et al 

lyalentjnTet al. 



V03 
Jl; K04b 



DP02; 



Di Paola et al. 



1I2OOOI): K OO: iKrisciunas et al.l 

I 2003h; K03: Krisciunas et al.l ||200 ; K04a: IKrisciunas et al, 
Krisciuna s et al. I ll2004d): K07: " ^CrkciunasetaL I ||2007^: ER06 : 



lElias-Rosa et al.l 1I2OO 6'); ER07: lElias-Rosa et all ll2007h : Pa07b: IPastorello et al.l ^OOJ^; 
StOT: IStanishev et al.l li2007l) ; P08: IPienata et al.l ll2008[l . 

•^The fast-declining, unusual SN 2005bl, SN 2005hk, and SN 2005ke were excluded from 
the template generation. When fit with the If-band template they are clearly rejected as 



these dromedary SNe la compared with the H-band template. 

''Because we had data from PAIRITEL on these SNe la we did not take advantage of 
the photom etry of SN 2005cf provided in IPastorello et al.l bOOTd) or that of SN 2006X 
provided in I Wang et"al] 1I2OO8I) . 
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Table 2. SN la Ji/i^^Templates 



Phase J Magnitude Uncertainty H Magnitude Uncertainty Ks Magnitude Uncertainty 



-10.00 


0.59 


0.26 


0.45 


0.23 


0.64 


0.26 


-9.00 


0.36 


0.30 


0.26 


0.26 


0.51 


0.33 


-8.00 


0.24 


0.25 


0.16 


0.23 


0.33 


0.27 


-7.00 


0.10 


0.20 


0.04 


0.21 


0.20 


0.21 


-6.00 


0.00 


0.13 


-0.03 


0.18 


0.10 


0.09 


-5.00 


-0.07 


0.09 


-0.08 


0.14 


0.03 


0.10 


-4.00 


-0.11 


0.07 


-0.10 


0.13 


-0.02 


0.11 


-3.00 


-0.11 


0.07 


-0.09 


0.13 


-0.04 


0.10 


-2.00 


-0.09 


0.09 


-0.07 


0.13 


-0.04 


0.10 


-1.00 


-0.06 


0.11 


-0.04 


0.13 


-0.03 


0.10 


0.00 


0.00 


0.11 


0.00 


0.12 


0.00 


0.13 


1.00 


0.07 


0.13 


0.04 


0.11 


0.04 


0.14 


2.00 


0.16 


0.15 


0.07 


0.10 


0.08 


0.14 


3.00 


0.26 


0.19 


0.10 


0.12 


0.12 


0.13 


4.00 


0.37 


0.24 


0.13 


0.13 


0.15 


0.13 


5.00 


0.50 


0.26 


0.15 


0.12 


0.19 


0.13 


6.00 


0.64 


0.25 


0.18 


0.12 


0.22 


0.12 


7.00 


0.78 


0.26 


0.20 


0.14 


0.25 


0.14 


8.00 


0.92 


0.27 


0.23 


0.14 


0.28 


0.14 


9.00 


1.05 


0.26 


0.25 


0.13 


0.31 


0.17 


10.00 


1.17 


0.26 


0.26 


0.13 


0.33 


0.18 


11.00 


1.28 


0.24 


0.25 


0.13 


0.35 


0.17 


12.00 


1.37 


0.22 


0.24 


0.11 


0.35 


0.18 


13.00 


1.43 


0.19 


0.22 


0.12 


0.34 


0.17 


14.00 


1.47 


0.16 


0.19 


0.13 


0.33 


0.15 


15.00 


1.48 


0.15 


0.15 


0.13 


0.31 


0.14 


16.00 


1.47 


0.16 


0.12 


0.11 


0.30 


0.17 


17.00 


1.44 


0.19 


0.09 


0.11 


0.29 


0.16 


18.00 


1.40 


0.20 


0.06 


0.11 


0.28 


0.17 


19.00 


1.36 


0.23 


0.03 


0.13 


0.26 


0.18 


20.00 


1.32 


0.26 


0.00 


0.13 


0.24 


0.18 


21.00 


1.26 


0.26 


-0.03 


0.12 


0.22 


0.16 


22.00 


1.21 


0.23 


-0.06 


0.11 


0.21 


0.17 


23.00 


1.17 


0.27 


-0.07 


0.12 


0.22 


0.17 


24.00 


1.11 


0.29 


-0.07 


0.12 


0.24 


0.17 


25.00 


1.06 


0.29 


-0.06 


0.14 


0.26 


0.19 


26.00 


1.02 


0.28 


-0.04 


0.16 


0.30 


0.19 


27.00 


0.98 


0.26 


0.00 


0.22 


0.33 


0.21 


28.00 


0.96 


0.26 


0.03 


0.24 


0.35 


0.21 


29.00 


0.96 


0.25 


0.06 


0.28 


0.36 


0.21 


30.00 


0.97 


0.24 


0.10 


0.29 


0.38 


0.22 


31.00 


1.00 


0.25 


0.13 


0.25 


0.40 


0.25 


32.00 


1.04 


0.25 


0.17 


0.25 


0.43 


0.25 


33.00 


1.10 


0.25 


0.21 


0.25 


0.47 


0.25 


34.00 


1.16 


0.25 


0.26 


0.25 


0.52 


0.25 
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Table 2 — Continued 



Phase 


J Magnitude 


Uncertainty 


H Magnitude 


Uncertainty K. 


j Magnitude 


Uncertainty 


35.00 


1.23 


0.25 


0.30 


0.25 


0.57 




0.25 


36.00 


1.30 


0.25 


0.35 


0.25 


0.62 




0.25 


37.00 


1.37 


0.25 


0.40 


0.25 


0.67 




0.25 


38.00 


1.44 


0.25 


0.45 


0.25 


0.72 




0.25 


39.00 


1.51 


0.25 


0.49 


0.25 


0.77 




0.25 


40.00 


1.58 


0.25 


0.54 


0.25 


0.82 




0.25 


50.00 


2.27 


0.25 


1.01 


0.25 


1.31 




0.25 


60.00 


2.97 


0.25 


1.49 


0.25 


1.80 




0.25 


70.00 


3.67 


0.25 


1.96 


0.25 


2.30 




0.25 


80.00 


4.37 


0.25 


2.44 


0.25 


2.79 




0.25 


Note. 


— The templates are linear 


fits from +35 


days to +80 days 


with (J, H, 


Ks) s 


lopes of 



(0.0625, 0.0470, 0.425) mag/day. 



Table 3. Apparent Standard if-Band Magnitudes of Type la Supernovae 



SN Name 


C^CMB+Virgo'' CZeii 






X^/DoF DoF'l Ay 








km km 


[mag] [MJD] 


[mag] [mag] 


[mag] 


[mag] 





SN 1998bu 


709 


20 


0.463 


50952.4 


11.76 


0.04 


0.42 


12 


1.12 


29.76 


0.46 


SN 1999cl 


957 


86 


0.392 


51342.2 


12.76 


0.06 


1.42 


3 


3.49 


30.77 


0.40 


SN 1999cp 


2909 


14 


0.112 


51363.2 


14.89 


0.11 


0.20 


1 


0.08 


32.90 


0.16 


SN 1999ee 


3296 


15 


0.099 


51469.3 


15.11 


0.03 


0.14 


13 


0.96 


33.12 


0.11 


SN 1999ek 


5191 


10 


0.063 


51481.8 


16.35 


0.04 


0.30 


12 


0.54 


34.36 


0.07 


SN 1999gp 


8113 


18 


0.040 


51550.1 


17.07 


0.10 


0.24 


2 


0.10 


35.08 


0.11 


SN 2000E 


1803 


19 


0.182 


51577.2 


13.84 


0.05 


0.51 


8 


0.66 


31.84 


0.19 


SN 2000bh 


6765 


21 


0.049 


51636.0 


16.85 


0.04 


0.78 


7 


0.16 


34.86 


0.06 


SN 2000bk 


7976 


20 


0.041 


51647.7 


17.35 


0.03 


0.50 


12 


0.30 


35.36 


0.05 


SN 2000ca 


6989 


62 


0.050 


51666.2 


16.78 


0.04 


0.74 


9 


0.03 


34.79 


0.06 


SN 2000ce 


5097 


15 


0.064 


51667.1 


16.13 


0.07 


0.36 


2 


1.87 


34.14 


0.09 


SN 2001ba 


8718 


22 


0.038 


52034.2 


17.52 


0.05 


0.51 


8 


0.05 


35.53 


0.06 


SN 2001bt 


4220 


13 


0.077 


52062.9 


15.84 


0.05 


0.42 


9 


0.69 


33.84 


0.09 


SN 2001cn 


4454 


250 


0.142 


52071.0 


15.82 


0.05 


0.25 


6 


0.52 


33.83 


0.15 


SN 2001CZ 


4506 


20 


0.073 


52103.4 


15.84 


0.07 


0.05 


6 


0.30 


33.85 


0.10 


SN 2001el 


978 


10 


0.334 


52182.1 


13.00 


0.05 


0.38 


8 


0.79 


31.01 


0.34 


SN 2002bo 


1696 


20 


0.194 


52356.0 


13.92 


0.05 


0.15 


6 


1.53 


31.93 


0.20 


SN 2002CV 


1696 


20 


0.194 


52417.7 


13.87 


0.06 


0.43 


7 


8.74 


31.88 


0.20 


SN 2002dj 


2880 


22 


0.114 


52450.6 


14.86 


0.04 


0.20 


10 


0.41 


32.87 


0.12 


SN 2003cg 


1340 


24 


0.246 


52729.1 


13.56 


0.07 


0.03 


4 


4.20 


31.57 


0.25 


SN 2003du 


2206 


14 


0.148 


52766.2 


14.62 


0.05 


0.60 


4 


0.04 


32.63 


0.16 


SN 2004S 


2607 


16 


0.126 


53038.7 


14.95 


0.05 


0.24 


5 


0.40 


32.96 


0.13 


SN 2004eo 


4859 


17 


0.067 


53278.7 


15.91 


0.07 


0.30 


4 


0.01 


33.92 


0.10 


SN 2005ao 


11828 


126 


0.036 


53442.0 


18.10 


0.09 


0.08 


1 


0.10 


36.11 


0.09 


SN 2005bl 


7360 


29 


0.045 


53482.1 


17.44 


0.05 


3.65 


9 


0.22 


35.45 


0.07 


SN 2005cf 


2018 


11 


0.162 


53533.6 


14.00 


0.04 


0.14 


7 


0.29 


32.01 


0.17 


SN 2005ch 


8094 


1499 


0.404 


53536.0 


17.26 


0.06 


0.73 


5 




35.27 


0.41 


SN 2005el 


4349 


8 


0.075 


53646.1 


15.79 


0.03 


0.78 


14 


0.03 


33.80 


0.08 


SN 2005eq 


8535 


25 


0.039 


53653.9 


17.34 


0.04 


2.11 


11 


0.11 


35.35 


0.06 


SN 2005eu 


10503 


14 


0.031 


53659.8 


17.38 


0.05 


0.87 


11 


0.08 


35.39 


0.06 


SN 2005hk 


3855 


21 


0.085 


53685.7 


15.68 


0.04 


13.09 


8 


1.19 


33.69 


0.10 


SN 2005iq 


10102 


40 


0.033 


53687.1 


17.73 


0.15 


0.06 


2 


0.05 


35.74 


0.15 


SN 2005ke 


1230 


11 


0.266 


53699.3 


14.05 


0.06 


8.32 


5 


0.34 


32.06 


0.27 



Table 3 — Continued 



oiN IN aiiie 


a 

C2;CMB+Virgo 


CZcTT 


fj ^ {VCL ) 


f 7-> 


^iBmax 


^(.-HBmaxJ 






AV 








km s ^ 


km s 


[mag] 


[IVi J UJ 


[mag] 


[mag] 


[mag] 






[mag] 




SN 2005na 


7826 


26 


0.042 


53740.5 


17.23 


0.10 


0.28 


5 


0.08 


35.24 


0.11 


SN 2006D 


2560 


18 


0.128 


53756.7 


14.67 


0.03 


0.33 


12 


0.11 


32.68 


0.13 


SN 2006N 


4468 


27 


0.074 


53760.6 


15.94 


0.12 


0.23 


4 


0.12 


33.95 


0.15 


SN 2006X 


1091 


20 


0.301 


53785.5 


12.79 


0.03 


0.24 


17 


3.83 


30.80 


0.30 


SN 2006ac 


7123 


17 


0.046 


53781.2 


16.88 


0.08 


0.23 


13 


0.08 


34.89 


0.09 


SN 2006ax 


4955 


20 


0.066 


53826.7 


16.36 


0.05 


0.38 


9 


0.14 


34.37 


0.09 


SN 2006cp 


6816 


14 


0.048 


53896.7 


16.83 


0.07 


0.63 


4 


0.75 


34.84 


0.09 


SN 2006gr 


10547 


22 


0.031 


54012.2 


18.12 


0.15 


1.16 


3 


0.19 


36.12 


0.15 


SN 20061e 


5403 


12 


0.060 


54047.2 


16.64 


0.05 


0.38 


16 


0.09 


34.65 


0.08 


SN 20061f 


4048 


10 


0.081 


54044.8 


15.81 


0.05 


0.18 


21 


0.26 


33.82 


0.09 


SN 2007cq 


7501 


50 


0.046 


54272.5 


16.97 


0.14 


1.26 


3 




34.98 


0.15 



''Redshift of host galaxies as corrected to CMB frame with the additional correction of Virgo member NGC 4501 (host of SN 1999cl) ^ 
from its observed CMB-corrected recession velocity (2281 km/s) to the mean recession velocity of Virgo (957 km/s). cn 

''A peculiar velocity of 150 km/s and the individual redshift measurement uncertainty of column 2 was converted into an equivalent 
distance modulus uncertainty using Eq. [T] 

'^An H-band extinction correction of Ah = O.OdAy has been applied to the magnitudes obtained by fitting to the template. See ^ 
for details. 

''DoF: The number of degrees of freedom is the number of //-band data points minus 1 for the overall offset fit parameter, i/emax- 
^Assuming Hq = 72 km s'^ Mpc-^and thus Mh = -18.07 ± 0.03 mag. 

'Quadrature sum of peculiar velocity distance modulus uncertainty (column 4) and fit uncertainty (column 7). 
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Table 4. JHK Absolute Magnitudes of SNe la at tsmax 



Filter PAIRITEL sample Global sample 

M RMS M RMS 



J -18.29 ±0.09 0.33 -18.27 ±0.07 0.29 

H -18.08 ±0.04 0.15 -18.07 ±0.03 0.16 

K -18.32 ±0.11 0.26 -18.30 ±0.11 0.29 



Note. — The J-band template shows significant varia- 
tion as a function of optical light-curve width. The num- 
bers derived above are for the average template as shown 
in Fig. [31 The K-band data from PAIRITEL is only useful 
for 9 new SNe la and does n ot sign ificantly improve on the 
results of iKrisciunas et al.l (l2004al ). 



Table 5. Light Curve Table Stub 



SN MJD Passband Fluxss^ Flux erras 



2005ao 53440.4157 H 

2005ao 53456.4300 H 

2005ao 53464.3872 H 

2005ao 53466.3879 H 

2005ao 53467.3847 H 



600.2491 3.5737 

437.9836 3.6405 

339.2386 128.0434 

345.1047 90.6351 

561.7854 3.7594 



^Fluxes are expressed normalized to a zeropoint of 
25. I.e., the cahbrated 2MASS magnitude is mag= 
-2.5 logio(flux) + 25 

Note. — This is a representative stub of the light-curve 
table. JHKs light curves for all 21 PAIRITEL SNe la pre- 
sented in this paper are available in the online electronic 
version. 



